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Abstract 

Controlling the type and density of charge carriers by doping is the key step for developing 



c/3 

1^ ■ graphene electronics. However, direct doping of graphene is rather challenge. Based on first- 



principles calculations, a concept of overcoming doping difficulty in graphene via substrate is re- 
ported. We find that doping could be strongly enhanced in epitaxial graphene grown on silicon 



'^ ■ carbide substrate. Compared to free-standing graphene, the formation energies of the dopants can 

O . decrease by as much as 8 eV. The type and density of the charge carriers of epitaxial graphene 



layer can be effectively manipulated by suitable dopants and surface passivation. More impor- 



^ ■ tantly, contrasting to the direct doping of graphene, the charge carriers in epitaxial graphene layer 

VO . are weakly scattered by dopants due to the spatial separation between dopants and conducting 

\0 ■ 

'^ I channel. Finally, we show that similar idea can also be used to control magnetic properties, e.g., 

\>- 

i^ ' induces a half-metallic state in the epitaxial graphene without magnetic impurity doping. 



PACS numbers: 



The conventional silicon based microelectronics is expected to encounter fundamental 
limitations at nanoscale. According to the semiconductor industry road map, novel materials 
that could complement or substitute silicon are needed. It is shown that graphene is suitable 
for coherent nanoscale electronics applications due to its unique electronic properties and 
extremely high carrier mobility [l|-l3|. Controlling the type and density of charge carriers 
by doping is at the heart of graphene electronics. However, the development of reliable 
chemical doping methods in graphene is still a real challenge. Direct doping of graphene 
with substitutional dopants is rather difficult, because the strong covalent C-C a bonds make 
graphene as one of the strongest materials in the world l|. Until now, doping of graphene has 



Q,Q, 



mainly been achieved through reactive molecular or atomic adsorbates in experiments 
which is difficult to implement under current device technology. Moreover, direct doping 









usually suppresses the high carrier mobility of graphene evidently |l|-|3|]. Thus, new proposals 
for overcoming the doping difficulty in graphene is urgently needed. 

Comparing to mechanical exfoliation of graphene from graphite, large area epitaxial 
growth of graphene layers on silicon carbide (SiC) surface shows huge application potentials 



3,E|. 



In 



and has apparent technological advantages over the mechanical exfoliation method 
this letter, we demonstrate that substitutional doping could be strongly enhanced in epitax- 
ial graphene grown on SiC substrate comparing to that of free-standing graphene (FSG). 
The formation energies of dopants can reduce by as much as ~ 8 eV. The type and density of 
carriers in epitaxial graphene layer could be effectively manipulated by suitable dopants and 
surface passivation. Contrasting to direct doping of FSG, the carriers of epitaxial graphene 
are weakly scattered by dopants due to the spatial separation between scattering centers 
(dopants) and conducting channel (epitaxial graphene). Thus, the intrinsic high carrier mo- 
bility of graphene could be effectively maintained even during doping process. Finally, we 
show that the same idea also can be used to control spin properties of epitaxial graphene by 
buffer layer defects. We find that the reconstructed defects in buffer layer can break the spin 
symmetry of epitaxial graphene and half-met allicity can be induced in the system without 
magnetic impurity doping. 

All the calculations are performed using the method based on density functional theory 
in the local spin density approximation (LSDA) as implemented in the VASP code|l8|, which 
is known to describe the structural and electronic properties of graphite quite welljol, Il0 |. 



The LSDA results are further checked by the generalized gradient approximation and the 
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FIG. 1: (color online) The side view of optimized structure of graphene layers on SiC(OOOl) surface. 
The top epitaxial graphene layer and the rippled interface buffer layer are highlighted in green/light 
gray and red/dark gray, respectively. Only half of the SiC substrate is plotted here. 

results are consistent with each other. The electron-ion interaction is described by the 
projector augmented wave method, and the energy cutoff is set to 400 eV. Structural opti- 
mization is carried out on all systems until the residual forces are converged to 0.02 eV/A. 
A v3 X \/3R30° SiC substrate cell is used to accommodate a 2 x 2 epitaxial graphene 
cell, approximating to the larger reconstruction with a 6v3 x 6v3i?30° periodicity usually 
observed in experiments. This approximation has been proved to well describe the elec- 
tronic properties of epitaxial graphene on SiC|9|, |ll|. More importantly, the a/3 x ^/3R30° 
reconstruction is also been observed in experiments during epitaxial growth of graphene on 
both (C-terminated) SiC(OOOl) surface [6|,[l2(] and (Si-terminated) SiC(OOOl) surface IsNlSj. 
To address the doping effects in graphene-SiC system, a 2v^ x 2-\/3/?30° SiC substrate cell, 
which accommodates a 4 x 4 graphene cell, is constructed. The SiC substrate is modeled 
by a slab contains 4 SiC bilayers with H passivation on the second surface of the slab. Sev- 
eral test calculations on 6 bilayers and 8 bilayers substates give essentially the same results. 
A gamma centered 6x6x1 k-point sampling is used for the Brillouin-zone integration, 
including the Dirac point. In the following, although our discussion is focused on the results 
of graphene on SiC(OOOl) surface, the results for graphene on SiC(OOOl) surface will also be 
mentioned. 

The optimized structure of graphene layers on SiC(OOOl) is shown in Fig. 1. 2/3 of 
the C atoms in the top SiC surface make covalent bonds with the C atoms of the interface 
buffer layer (the first graphene layer during vacuum graphitization on SiC), with a C-C 

y ch anges the Dirac-like band structures of the buffer 
ll| . A similar buffer layer structure is also observed 
on SiC(OOOl) surface. The second graphene layer on the top of buffer layer is planar and the 
distance between the rippled buffer layer and top graphene layer varies between 3.34 - 3.71 



bond length ~ 1.62 A, which drastical 

n 

layer, agreeing with previous results [9 



A. The interface buffer layer is stretched by ~ 8% as a result of the lattice mismatch between 
2-\/3 X 2\/3R30° SiC surface and graphene, agreeing with previous predictions |9l-lll|. 



It is expected that the carriers in graphene could be manipulated by p/n-type dopants, 



like in conventional semiconductors |l6l|. Al, B, P, and N are considered as substitutional 
dopants, as they are widely used in C- and Si- based materials. All the possible substitutional 
positions in graphene-SiC system are checked and the formation energies of dopants at 
different positions are calculated to search the most stable configurations. The formation 
energy (Ef) of a substitutional dopant in graphene-SiC system is defined as: 

Ef = Edoped — Ehost + IJ-C/Si " IJ-dopant (1) 

where Eloped and Ehost are the total energies of doped and undoped system, respectively. 
fic/Si and Udopant (dopant = Al, B, P, N) are the chemical potentials of C (or Si) and dopant, 
respectively. Since epitaxial graphene should grow in C-rich environment, ^c is calculated 
assuming epitaxial graphene is stable and usi is calculated so that SiC is stable. It should 
be noticed that the relative stability between different doping positions does not depend on 
the particular choice of Udopant- The calculated formation energy differences are shown in 
Fig. 2. Interestingly, Al prefers to substitute Si atom of the SiC substrate or substitute 
C atom of the interface buffer layer and the Ef of Al decreases dramatically by ~ 8.26 eV 
compared to FSG case. B and P atoms prefer to substitute C atoms of the buffer layer and 
the Ef of dopants decrease largely by 2.36 eV and 5.27 eV, respectively, compared to FSG 
cases. Besides doping in buffer layer, the Ef of B and P at SiC surface are also reduced by 
1.46 eV and 4.87 eV, respectively, compared to FSG cases. Contrary to the cases of Al, B, 
and P doping, N atom prefers to substitute C atom of top SiC substrate with a reduction of 
Ef by 3.52 eV, as shown in Fig. 2d. The similar phenomenon of substrate-enhanced doping 
is also found on SiC (0001) surface, with the Ef of dopants decreased by as much as 10.8 
eV. 

These quite different doping behaviors between different dopants as well as different 
layers in graphene-SiC system could be understood by considering the structural strain 
in the various epitaxial graphene-SiC layers, the different atomic covalent radiuses, and the 
electronegativity of dopants. Since the interface buffer layer is under large tensible strain, Al, 
B, and P atoms could be easier to substitute C atoms in this layer due to the larger covalent 
radiuses of Al (1.21 A), B (0.84 A), P (1.07 A) than s/ C (0.73 A). The rippled (deformed) 
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FIG. 2: (color online) The calculated formation energy differences of dopants (a) Al, (b) B, (c) 
P, and (d) N in graphene-SiC system as a function of different layers. The formation energies of 
dopants in free-standing graphene (FSG) are also plotted in these figures for comparison. 

structure reduces the doping difficulty in buffer layer to some extent. Furthermore, because 
the electronegativity of C (2.55) is larger than Si (1.90) in SiC, Si behaves as positively 
charged cation and C acts as negatively charged anion. Al atom has lower electronegativity 
(1.61) and a little larger covalent radiuses than Si (1.11 A), so it would prefer to substitute 
Si atom in the SiC matrix to increase the negative Coulomb interaction between Al and C. 
Differing from Al, B, and P, the structural tensile strain increases the difficulty of N doping 
in the graphene layer because N has a smaller covalent radius (0.71 A) than C. N atom 
has even higher electronegativity (3.04) than C, so it would prefer to substitute C atom 
in the SiC matrix to increase the negative Coulomb interaction between Si and N, which 
is in agreement with the calculated trend of Ef in Fig. 2d. More importantly, the large 
reduction of Ej for all the dopants at strain-free SiC substrate (1.46 ~ 8.26 eV) indicates 
that, not only in the 2v3 x 2v^-R30° structure, the strongly substrate-enhanced doping can 
also exist in 6-\/3 x 6\/3R30° reconstruction because the strain-free substrate is independent 
of reconstruction in the graphene or buffer layers. The above results and analysis show that 
overcoming the doping difficulty in graphene could be generally achieved through substrates, 
and it can also be applied to other layered or thin-film systems. Moreover, it is quite 



encouraging to see that highly strained graphene (vm 



substrates in experiments through various schemes 17|, 



to 25%) have already been created on 
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FIG. 3: (color online) The electronic band structures of graphene-SiC system without doping (a) 
and doped by (b) B and (c) P in buffer layer, (d)-(f) are the same as (a)-(c) but the C atoms with 
dangling bonds at the SiC surface are passivated by H atoms. The blue/dark gray solid line and 
green/light gray dotted line represent the spin-up and spin-down states, respectively. The Fermi 
level is set at zero energy. 



After knowing the stable doping configurations, we turn to investigate how these dopants 
tune the carrier properties, i.e. carrier type and carrier density, in graphene-SiC system. 
Without doping, the Dirac-hke states in the top epitaxial graphene layer on SiC (0001) is 
preserved, as shown in Fig. 3a. The weakly dispersive interface states visible in Fig. 3a 
result from the surface dangling bonds of SiC substrate. In our unit cell, there are twelve 
C sp^ dangling orbitals on the top layer of the SiC substrate. Eight of them make covalent 
bonds with the buffer layer but the other four remain unsaturated. The localization of C sp^ 
dangling bonds favor the spin-polarization and thus split to four occupied states and empty 
states, respectively. Because the band alignment between the epitaxial graphene states and 
substrate surface states is type-I, therefore, there is no charge transfer between the epitaxial 
graphene and substrate surface. However, the following results show that the substrate 
surface states play an important role in affecting the carrier density of top graphene layer. 

Our total energy calculations indicate that dopants favor to locate at buffer layer or c 
to buffer layer in the substrate. To understand the carrier distribution, Bader analysis 
is used to estimate the charge variation in the dopant and various layers. We find that 



ose 

3 



after substitutional doping in buffer layer, B atom shows an onsite charge of 1.32 e, whereas 
in the case of Al doping, the value is substantially lower, 0.80 e. The variation between 
the Al and B cases could be understood by noticing that B is much more electronegative 
than Al. When the total hole carrier density is ~ 1.02 x 10^^ cm~^, the hole density in 
the buffer layer is ~ 5.77 x 10^^ cm^^ after B doping and ~ 9.98 x 10^'^ cm^^ after Al 
doping. Unexpectedly, the top epitaxial graphene is not effectively doped by holes in either 
cases, with hole density at ~ 5 x 10^^ cm~^. The origin of the small carrier density in the 
epitaxial graphene layer can be understood from the band structure of the graphene-SiC 
system. We see that there are localized interface states which strongly pin the Fermi level at 
about 0.18 eV below the Dirac point for ^type materials, as shown in Fig. 3b. (Al-doping 
case is similar and not shown here). The hole density of epitaxial graphene is thus limited 
to ~ 5 X 10^^ cm~^ for p-tjpe doped system. For the n-doped case (Fig. 3c), the empty 
interface states also play a similar role and pin the Fermi level at about 0.44 eV above the 
Dirac point, which is responsible to the limited electron density at ~ 1.99 x 10^^ cm~^ in 
the top epitaxial graphene layer after P doping. The above results strongly demonstrate 
that although B (Al) and P doping in buffer layer could induce electrons or holes in the top 
epitaxial graphene, the interface states play as a " subthreshold valve" to restrict the carrier 
density. Similar doping behaviors are found when B (Al) and P dope at SiC surface. As 
we mentioned, N atom prefers to substitute the C atom with dangling bond at the top SiC 
layer. Contrary to P doping, the epitaxial graphene layer could not be doped with electrons 
because the N atom is strongly electronegative so the spin-split states of N atom are both 
occupied after doping. 

An obvious route to improve the electron or hole density in epitaxial graphene is by 
eliminating the influence of the interface states. Since these interface states originate from 
the dangling orbitals of C atoms at the top of SiC substrate, naturally, we may eliminate 
them by surface passivation. Here, we take H as an example to demonstrate the surface 



passivation effect on modulating the charge density of epitaxial graphene. In facL 



i in- 
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m- 



tercalation in epitaxial graphene has already been achieved in several experiments 
The main effect of H adsorption is the disappearance of the interface states, as shown in 
Fig. 3d. Without these interface states, we expect more electrons or holes will transfer from 
buffer layer to epitaxial graphene layer. The band structure calculations confirm our idea, 
as shown in Fig. 3e and 3f for B and P doping, respectively. After surface H passivation, the 
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FIG. 4: (color online) (a) The top view of optimized structure of graphene-SiC system with recon- 
structed vacancy in the buffer layer (small red/dark gray atoms), (b) The electronic band structure 
of (a). The blue/dark gray solid line and green/light gray dotted line represent the spin- up and 
spin-down states, respectively. The Fermi level is set at zero energy. 



hole density of B-doped system increases four times up to 2.48 x 10^'^ cm^^, which pushes 
the the Fermi level to 0.42 eV below the Dirac point. Similarly, the electron density of P- 
doping system increases twice to 3.58 x 10^^ cm~^. Moreover, a small band gap (~ 0.07 eV) 
appears in P-doping system, as shown in Fig. 3f. These results strongly demonstrate that 
not only we can control the type of carriers in epitaxial graphene via doping of the buffer 
layer, but also the carrier density could be controlled by substrate passivation. More notice- 
ably, the dopants prefer to stay in the buffer layer while the carriers are concentrated in the 
top epitaxial graphene. Thus, the carriers in epitaxial graphene layer could be only weakly 
scattered by dopants due to the spatial separation between scattering centers (dopants) and 
conducting channel (top graphene layer). This is in the same spirit as in modulation doping 
in heterostructures such as GaAlAs/GaAs[2.2|, and we suggest that it is an efficient way in 
maintaining the high carrier mobility of graphene during doping process. 

The manipulation of buffer layer not only can be used to control the charge carriers, 
but also could be used to tune the spins of the system. Here we demonstrate that the 
reconstructed vacancy in buffer layer has unexpected effect on the spin properties of epi- 
taxial graphene. Injection of high spin-polarized current in graphene is the current subject 
of intense investigation efforts J23|. Ideally, 100% spin-polarized current could be induced 
by half-metallic materials. Although spontaneous vacancy in FSG may induce local spin- 
polarization, it is rare due to its high formation energy of ~ 8 eV. In graphene-SiC system, 
we find that vacancy prefers to be formed in buffer layer and undergoes reconstruction, as 
shown in Fig. 4a. A five-thirteen ring is formed with a Ef of 3.2 eV, which is ~ 5 eV lower 
than in FSG. The deformed buffer layer as well as the defect reconstruction are responsible 



for the lower Ej of vacancy in buffer layer. In experiments, considerable amounts of defects 
have been observed in the interface layer of epitaxial graphene system [ij, |2J], especially the 



existence of reconstructed hexagon-pentagon-heptagon 2J]. The experimental observations 
indicate that the Ef of defect in buffer layer is much lower than in FSG, agreeing with our re- 
sults. It is interesting to see that the reconstructed vacancy not only breaks the honeycomb 
symmetry, but also breaks the spin symmetry of graphene, which induces a half-metallic 
state in the top epitaxial graphene layer, as shown in Fig. 4b. There is an apparent gap 
(~ 0.142 eV) in the spin-up state of the epitaxial graphene and a negligible band gap in 
spin-down state. The localized states near the band edge in spin-up band structure are 
induced by the C atoms which are bonded to the SiC surface in the five-thirteen ring. It 
was theoretically predicted that half-metallicity could be induced in graphene by an exter- 



nal electric field 25(] or transition metal doping 26| . However, the required electric field is 



too strong and heavy transition metal elements also often act as poison agents in biological 
systems. Our results here strongly indicate that SiC-graphene system with reconstructed 
defects in the buffer layer could have strong potentials for spintronics. 

In conclusion, a concept of overcoming doping difficulty in graphene via substrate is re- 
ported. We show that substitutional doping could be strongly enhanced in epitaxial graphene 
on SiC substrate. Compared to free-standing graphene, the formation energies of dopants 
decrease by as much as 8 eV. The type and density of carriers in epitaxial graphene could be 
effectively manipulated by suitable dopants and surface passivation. The carriers in epitaxial 
graphene layer are weakly scattered by dopants due to the spatial separation between the 
dopants and carriers. Finally, we show that the reconstructed vacancy in buffer layer could 
induce half-metallicity in epitaxial graphene without magnetic impurity doping. Generally, 
the effect of substrate-enhanced doping could exist in other substrates, and it can also be 
applied to other layered or thin-film systems. 
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No. DE-AC36-08GO28308. The work at Fudan was partially supported by the National 
Science Foundation of China. 
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